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/* amfm_8bit.c file */ 

/* 8 bits/pixel am/fm halftoning algorithm (with partial doting) */ 
/* Image length and width are assumed to be multiple of 8 */ 
/* One row serpentine TDED with suppressing each other dot */ 
/* The output is a tiff file containing 8bit pwm codes */ 

<stdio . h> 
<stdlib . h> 
<math . h> 
< t ime . h> 
"coef.h" 
"tiff . h" 
"allocate .h" 

void amfm (unsigned int, unsigned int, unsigned char *, unsigned char *,\ 

short *, unsigned char **,TDEDPARA *, short *, short *); 
void amfm_ed_8bits (unsigned int, unsigned int, unsigned char **, unsigned char**,\ 

unsigned char **,TDEDPARA *, short *, short *); 
int main (int argc, char ** argv) 

{ 

int i , j ; 

unsigned int height , width ; 
FILE * fp; 

struct TIFF_img input__img, output_img, mid; 
time_t first, second; 

TDEDPARA *tdedpara = &TDEDcoef f [0 ] ; 

short *dotdensityLUT = &OptDensityLUT [0] ; 

short *dotsizeLUT = &OptSizeLUT [0 ] ; 

if(argc<3) { 

printf ( "usage : %s input_img . tif output_img . tif \n" , argv [0] ) ; 
return 1; 

} 

/* read the input image */ 

if ( {fp = fopen (argv [1] , " rb" ) ) ==NULL) { 



y=; #include 

Q #include 

= #include 

fl #include 

ffi #include 

5l #include 
#include 




printf("can not open file %s l\n" , argv [1] ) ; 
exit (2) ; 

} 

if (read_TIFF (fp, &:input_img) ) { 

printf ( "error reading input file\n"); 
exit (3) ; 

} 

f close ( f p) ; 

if ( (fp=fopen ( "dbshalf . tif " , "rb" ) ) ==NULL) { 

fprintf (stderr, "can not open file: dbshalf . tif\n" ) ; 
exit (1) ; 

} 

if (read_TIFF(fp, &mid) ) 
{ 

fprintf (stderr , "error reading file\n"); 
exit (1) ; 

} 

f close (f p) ; 

/* Set variable to do timing of algorithm */ 
first = time (NULL); 

/* Modify image width to make sure each strip is multiple of 8 */ 
width = floor ( input_img . width/ 8) *8 ; 
height = floor { input_img . height/8) *8 ; 

/* Allocate memory for entire fm output image. */ 
get_TIFF( &:Output_img , height, width, *g' ); 

amf m_ed_8bits (height , width , input_img . mono , output_img . mono , \ 
mid . mono , tdedpara, dotdensityLUT, dotsizeLUT) ; 

/* show the run time */ 
second = time (NULL) ; 

fprintf (stdout ," \nFinished AM/FM and writing results . \n" ) ; 

f printf (stdout , "Cum - run time: %f sec . \n" , diff time (second, first )) ; 

/* write PWM codes image */ 

if( (fp = fopen(argv [2] , "wb") ) ==NULL) { 

printf ("cannot open file %s\n" , argv[2]); 

exit (4) ; 

} 

if (write_TIFF ( f p , &output_img) ) { 

printf ("\nError writing TIFF file %s\n" , argv[2]); 
return 1 ; 

} 

f close (fp) ; 

/* free the space */ 

f ree_TIFF (& (output_img) ) ; 

f ree_TIFF (& (input_img) ) ; 

f ree_TIFF (& (mid) ) ; 

f flush (stdout) ; 

return 0 ; 

} 



void amf m_ed_8bits ( 

unsigned int height, /* 
unsigned int width, /* 
unsigned char ** contone_img , 



unsigned char ** tok:en_img, 
unsigned char ** dbs_screen, 
part */ 

TDEDPARA *tdedpara, 
short * dotdens i t yLUT , 
short *dotsizeLUT) 



Input image height */ 
Input image width */ 

/* Input image [height] [width] */ 
Output token image [height] [width] */ 

/* DBS screen used in thresholding of fm 

Tone -dependent error diffusion parameters */ 
Optimal dot density curve */ 
Optimal dot size curve */ 



{ 



short *fm_err; 
unsigned int i/j; 

/* initialize first row of fm error buffer */ 
srand(l); /* fix the seed */ 

fm_err = (short* ) malloc (sizeof ( short ) * (width+2)); 
for(j = 0; j<width+2; j++) 

fm_err[j] = (rand ( ) %128 -64 ) ; /* initialization */ 

/* Process the input image with 2 rows each time */ 
for(i=0; i<height; i+=2) { 

if{(i%600) == 0) printf ( "amfm_ed: starting row %d\n", i) ; 

amf m (width, i , contone_img [i ] , token_img [i] , f m_err , dbs_screen, \ 
tdedpara, dotdens ityLUT, dot si zeLUT) ; 

} 

free (fm_err) ; 
return; 



} 



/* This subroutine only processes 2 rows */ 
/* Assume width of image is multiple of 8 */ 
void amf m ( 

unsigned int width, /* Input image width */ 
unsigned int i, /* ith row */ 

unsigned char *img_in, /* ith row of input image array */ 

unsigned char *img_out, /* ith row of output image array */ 

short *fm_err, /* FM error buffer */ 

unsigned char ** dbs_screen, /* dbs_screen [SCREENHEIGHT] [SCREENWIDTH] */ 

TDEDPARA *tdedpara, /* Tone-dependent error diffusion parameters */ 
short *dotdensityLUT, /* Optimal dot density curve */ 
short *dotsizeLUT) /* Optimal dot size curve */ 



short f m_tmp , thresholding ; 
short *f m_err_ptr , *tded_ptr ; 
short pixela, pixelb, output; 
int j ; 

unsigned char *img_in ptr, * img_out_pt r , *dbs_pat_rowptr ; 

short dotdensity, mod_input , error; 

short Wl, W2, W3, W4, T2, DT, el, e2 , e3 , e4 ; 

FILE *fp; 



/* serpentine even rows */ 
/* */ 

/* initial points */ 
f m_tmp = 0 ; 

fm_err_ptr = fm_err+l; 
img_in_ptr = img_in; 
img_out_ptr = img_out; 

dbs_pat_rowptr = dbs_screen [ {i++) %SCREENHEIGHT] ; /* Inverse dbs pattern */ 

/* Index through pixels in pairs */ 
for(j = 0; j<width; j-j+2 ) { 

/* First process FM (dot density) for left pixel in pixel pair. */ 

/* Get first pixel */ 
pixela = * (img_in_ptr++) ; 

/* Use look-up- table to get dot density */ 
dotdensity = dotdensityLUT [pixela] ; 

/* Compute look-up table entries for tone dependent error diffusion */ 
tded_ptr = (short*) (tdedpara + dotdensity) ; 



T2 = 


* (tded_j)tr++) 


DT = 


* (tded_j>tr++) 


Wl - 


* (tded_ptr++) 


W2 = 


* (tded_ptr++) 


W3 = 


* (tded_ptr++) 


W4 = 


*tded_j)tr; 



/* compute dotdensity modified by diffused error */ 
mod_input = dotdensity + *fm_err_ptr; 

/* Threshold modifed dotdensity */ 

thresholding = mod_input - (dbs j)at_rowptr [ j %SCREENWIDTH] * DT + T2) ; 
output = (thresholding > 0) ? 255 : 0; 

/* Compute weighted errors */ 
error = output - mod_input ; 
el = (Wl * error) >>8; 
e2 = (W2 * error) >>8; 
e3 = (W3 * error) >>8; 
/*e4^ = (W4 * error) >>8;*/ 
e4 = error - el - e2 -e3/ 

/* Diffuse error forward in 1-D error buffer */ 

* ( - - f m_err_ptr) -= e4 ; 

* (++fm_err_j)tr) = fm_tmp - e3 ; 

* (++fm_err_ptr) -= el; 

fm_tmp = -e2/ 

/* Now process FM (dot density) for right pixel in pixel pair. */ 
/* Use same TDED parameters as for left pixel. */ 

/* Get second pixel */ 
pixelb - * ( img_in_ptr++) ; 



/* Use look-up-table to get dot density */ 



dotdensity = dotdensityLUT [pixelb] ; 
mod_input = dotdensity + *fm_err_ptr; 

error = - mod_input; /* suppress dot firing at this pixel */ 

el = (Wl * error) >>8; 
e2 = (W2 * error) >>8; 
e3 := (W3 * error) >>8; 
/*e4 = (W4 * error) >>8;*/ 
e4 = error - el - e2 -e3; 

/* Using the tded weights of the left pixel */ 

* { - -fm_err__ptr) -= e4 ; 

* (++f m_err_ptr) = fm_tmp - e3 ; 

* (+ + fm_err j>tr) -= el; 

fm_tmp = -e2; 

/* Begin section on dot size rendering with partial doting */ 
if (output) { 

/* Left pixel */ 

* (img_out_ptr++) = (dotsizeLUT [pixela] >>1) +NEWRIGHT; 
/* Right pixel */ 

if (dotsizeLUT [pixela] & 1) /* Take care of quantization error */ 

* (img_out_ptr++) = ( (dotsizeLUT [pixelb] +1) >>1) + NEWLEFT; 
else 

* (img_out_ptr++) = (dotsizeLUT [pixelb] >>1) + NEWLEFT; 

} 

else { 

* (img_out_ptr++) = NEWRIGHT; 
* ( img_out_ptr++) = NEWLEFT; 

} 

} /* end of ith row */ 

/* */ 

/* serpentine odd rows */ 
/* */ 

f m_tmp = 0 ; 

/* Set fm error buffer pointer to the end of fm_err buffer */ 
fm_err_ptr = fm_err+ width - 1; /* offset by 1 */ 
img_in_ptr = img_in+width*2 -2 ; 
img_out_ptr = img_out+width*2 -1 ; 

*img_out_ptr = NEWRIGHT; /* Take care of the last pixel in odd row */ 
img_out_ptr -= 2; 

dbs_pat_rowptr = dbs_screen [i%SCREENHEIGHT] ; 

/* Index through pixels in pairs */ 
for(j = width-2; j> 0; j=j-2) { 

/* First process FM (dot density) for right pixel in pixel pair */ 

/* Get right pixel */ 

pixela = * ( img_in__ptr- - ) ; 

/* Use look-up- table to get dot density */ 
dotdensity = dotdensityLUT [pixela] ; 



/* Compute look-up table entries for tone dependent error diffusion */ 



tded_j)tr = (short*) (tdedpara + dotdensity) ; 



T2 = 


* (tded_ptr++) ; 


DT = 


* (tded_ptr+ + ) ; 


Wl = 


* (tded_ptr++) ; 


W2 = 


* (tded_ptr+ + ) ; 


W3 = 


* (tded_ptr++) ; 


W4 = 


*tded_j5tr; 



/* Compute dotdensity modified by diffused error */ 
mod_input = dotdensity + *fm_err_ptr; 

/* suppress this dot and compute the error */ 
error = - mod_input ; 

/* Compute weighted errors */ 

el = (Wl * error) >>8; 

e2 = (W2 * error )>>8; 

e3 = (W3 * error) >>8; 

/*e4 = ( (W4 * error) >>8) ; */ 

e4 = error - el - e2 -e3 ; 

/* duffuse error forward in 1-D error buffer */ 

* (++f m_err__ptr) -= e4 ; 

* ( - -f m_err_ptr ) = fm_tmp - e3 ; 

* ( - -f m__err_ptr ) -= el; 
fm_tmp = -e2; 

/* Now process FM (dot density) for Left pixel in a pair */ 

/* Get second pixel */ 
pixelb = * ( img_in_ptr- - ) ; 

/* Use look-up- table to get dot density */ 
dotdensity = dotdensityLUT [pixelb] ; 

mod_input = dotdensity + *fm_err_ptr; 

/* Threshold modi fed dotdensity */ 

thresholding = mod_input - (dbs_pat_rowptr [ ( j - 1) %SCREENWIDTH] * DT + T2); 
output = (thresholding > 0) ? 255 : 0; 

error = output - mod_input; 

el = (Wl * error) >>8; 
e2 = (W2 * error) >>8; 
e3 = (W3 * error) >>8; 
/*e4 = (W4 * error) >>8;*/ 
e4 = error - el - e2 -e3; 



* (++fm_err_j)tr) -= e4 ; 

* ( - -f m_err_ptr ) = fm_tmp - e3; 

* ( - -f m_err_ptr ) -= el; 
fm_tmp = -e2 ; 



/* Begin section on dot size rendering with partial doting */ 
if (output) { 



/ 



/* Left pixel */ 

* (img_out_ptr++) = (dotsizeLUT [pixelb] >>1) + NEWRIGHT; 
/* Right pixel */ 

if (dotsizeLUT [pixelb] & 1) /* Take care of quantization error */ 

*img_out_ptr = { (dotsizeLUT [pixela] +1) >>1) + NEWLEFT; 
else 

*img_out_ptr = (dotsizeLUT [pixela] >>1) + NEWLEFT; 

} 

else { 

* (img_out_ptr++) = NEWRIGHT; 
*img_out_ptr = NEWLEFT; 

} 

img_out_ptr -= 3; 

} 

* (++img_out_ptr) = NEWLEFT; /* Take care of the first column */ 
return; 

} 
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/* coef.h file */ 

#define SCREENHEIGHT 12 8 
#define SCREENWIDTH 12 8 

#define NEWRIGHT OxcO 
#define NEWLEFT 0x4 0 
#define NEWCENTER 0x00 



#define Fl 0x0007 
#define F2 0x0003 
#define F3 0x0005 
#define F4 0x0001 



/* Floyd-Steinberg Weights 7/16 in Q4 */ 

/* Floyd-Steinberg Weights 3/16 in Q4 */ 

/* Floyd-Steinberg Weights 5/16 in Q4 */ 

/* Floyd-Steinberg Weights 7/16 in Q4 */ 



typedef struct TDEDPARA 

{ 

short T2 
short DT 
short Wl 
short W2 
short W3 
short W4 
} TDEDPARA ; 



static TDEDPARA TDEDcoef f [12 9] = { 
{76, 0, 181, 0, 3, 72}, 
181, 0, 3, 72}, 
172, 1, 2, 81}, 
161, 14, 18, 63}, 
1, 37, 59}, 
6, 5, 96}, 
30, 0, 



{76, 
{79, 
{80, 
{82, 
{83, 
{83, 
{85, 
{86, 
{85, 
{86, 
{87, 



159, 
149, 
141, 
138, 
144 , 
129, 
123 , 
123, 



13, 
10, 
48, 
31, 
29, 



85}, 

105} , 

101} , 
76} , 

101} , 

101} , 



Of, 




115, 


28, 


5, 108}, 


Q Q 


U , 


138, 


19, 


18, 


81}, 


Q O 


r\ 


111, 


17, 


51, 


77} , 


Q Q 
O O , 


r\ 

0 / 


115, 


31, 


0, 110}, 


o / , 


rv 


120, 


16, 


16, 


104} 


Q Q 
O O , 


u , 


139, 


12, 


0, 105}, 




0 / 


122, 


19, 


17, 


98} , 


on 

yu / 


u , 


112, 


32 , 


0, 112}, 


Q T 

y ± , 




98, 34, 


20, 104}, 


yu , 


1 n 


123 , 


16, 26, 


91} 


y J # 


Q 

o , 


12 6 , 


1, 


74, 55}, 


O Q 

o y , 


T n 
X U , 


89 , 


26, 


71, 


70 } , 


o y , 


T n 


89, 


22 , 


43, 


102 } 


o y , 


1^ / 


91, 


21, 


34, 


110 } 


Q Q 
O O , 


xz , 


85, 


24, 


30, 


117 } 


o o 
c5 o , 


X4 , 


85, 


23, 


30, 


118 } 




A 

2 4 , 


113 , 


27, 13, 


103 




2, o , 


113 , 


33, 0, 


110} 


83 , 


2 6 , 


109, 


29, 9, 


109 } 


O A 


2 8 , 


106, 


21, 29, 


100 


o c 
ob , 


2 8 , 


103 , 


13, 56, 


84 } 


96 , 


2 , 


102 , 


16, 


57, 


81}, 


93 , 


6 , 


102 , 


25, 


28, 


101} 


91 , 


12 , 


102, 


24 


, 32, 


98 } 


96 , 


2 , 


103 , 


24, 


23, 


106 } 


94 , 


1 0 , 


99, 


17, 


62, 


78 } , 


95 t 


6 , 


110, 


12, 


110, 


24 } 


97 , 


4 , 


114 , 


12, 


112 , 


18 } 


y V , 


6 / 


114 , 


11. 


113 , 


18 } 


y D , 


o 

o , 


111, 


14, 


110, 


21 } 


O A 

y4 , 


X z , 


102 , 


17, 109 


, 2 8 


Q A 

y4 , 


o 


79, 32, 


108, 


37} , 


yb , 




74, 35, 


110, 


37} , 




2 , 


70, 35, 


111, 


40 } , 


y / , 


4 , 


68, 33, 


112, 


43) , 


y / , 


b , 


69, 28, 


112, 


47) , 


yo , 




70, 22, 


114, 


50) , 






68, 43, 


113 , 


) , 


1 U U , 


4 , 


68, 


22 , 


114, 


52} 


y y , 




71, 24, 


112 , 


49} , 


JL u ^ , 


z / 


70, 


23 , 


114, 


49} 


± u u , 


o , 


68, 


23, 


114, 


51} 


inn 

± L/ \J , 


Q 

O / 


66, 


22 , 


116, 


52 } 


1 no 


Q 


66, 


24 , 


116 , 


50} 




1 

X D , 


75, 


0, 


122, 


59}, 


y 3 , 


1 


63 , 


0, 


127 , 


66} , 


y 3 , 


Xb , 


56, 


0, 


130, 


70} , 




14 , 


56, 


0, 


132, 


68} , 


97 , 


16, 


59, 


0, 


132 , 


65}, 


97, 


16, 


60, 


0, 


133, 


63}, 


98, 


16, 


62 , 


0, 


133 , 


61}, 


95, 


26 , 


98, 


0, 


109, 


49} , 


97, 


20 , 


65, 


0, 


132 , 


59}, 


98, 


18 , 


61, 


0, 


132 , 


63}, 


99, 


18, 


63, 


0, 


131, 


62} , 


100, 


16, 58, 


0, 


133, 


65} 


100, 


16 


, 58, 


0, 


131, 


67} 



101 , 


16 , 


60 , 


0 , 


131, 


6 5), 


101 , 


16 , 


63 , 


0 , 


12 9, 


64 } , 


101 1 


16 , 


58 , 


0 , 


12 9, 


6 9}, 


102 , 


16 , 


7 1 , 


0 , 


12 3, 


62 } , 


103 , 


8 , 


68 , 


23 , 


114 , 


5 1 } , 


103 , 


8 , 


66 , 


22 , 


116, 


52 } , 


105 , 


6, 


68 , 


22 , 


115 , 


5 1 } , 


106 , 


4 , 


70 , 


22 , 


114 , 


5 0 } , 


108 , 


2 , 


69, 


23 , 


113 , 


5 1 } , 


105 , 


8 t 


68, 


22 , 


114 , 


52} , 


108 , 


6 , 


70, 


20, 


115, 


51} , 


106 , 


8 , 


69, 


27 , 


112, 


48} , 


109 , 


2 , 


65, 


35, 


112, 


44} , 


110 , 


4 , 


69, 


34, 


111, 


42} , 


110, 


6, 


72, 


35, 


110, 


39} , 


114 , 


0, 


73, 


34, 


111, 


38} , 


110, 


12, 


94 , 


21, 


108, 


33} 


111, 


12, 


102, 15, 110, 29 


116, 


6, 


114 , 


10, 


113, 


19} 


96, 


16, 


92, 


16, 


67, 81}, 


100, 


12, 


95, 


17, 


67, 


77} , 


101, 


12, 


97, 


19, 


67, 


73}, 



99, 4, 101, 20, 45, 90}, 

93, 4, 103, 25, 25, 103}, 

94, 8, 101, 25, 33, 97}, 
78, 24, 99, 26, 19, 112}, 

81, 26, 104, 22, 24, 106} 

82, 26, 102, 26, 25, 103} 
91, 26, 109, 14, 46, 87}, 

104, 10, 82, 0, 95, 79}, 
107, 8, 83, 0, 97, 76}, 

105, 8, 87, 2, 84, 83}, 
81, 14, 86, 27, 25, 118}, 
99, 12, 122, 0, 37, 97}, 

102, 10, 117, 0, 45, 94}, 

103, 10, 90, 21, 64, 81}, 
105, 12, 122, 4, 51, 79}, 
101, 12, 126, 9, 29, 92}, 

88, 12, 121, 25, 0, 110}, 

85, 12, 114, 25, 1, 116}, 

89, 10, 109, 23, 10, 114} 

86, 12, 112, 29, 1, 114}, 
89, 12, 119, 31, 0, 106}, 

94, 10, 123, 37, 1, 95}, 
93, 8, 117, 63, 1, 75}, 
99, 6, 118, 75, 9, 54}, 
97, 6, 120, 43, 3, 90}, 
111, 6, 121, 35, 32, 68}, 

95, 6, 116, 54, 0, 86}, 
107, 6, 125, 39, 15, 77}, 
93, 34, 137, 27, 19, 73}, 
85, 44, 139, 33, 16, 68}, 

87, 48, 146, 31, 23, 56}, 
87, 44, 148, 22, 10, 76}, 
93, 40, 152, 22, 11, 7l}, 
97, 44, 159, 4, 28, 65}, 
95, 42, 161, 25, 4, 66}, 



{103, 48, 176, 3, 44, 33}, 
{101, 56, 165, 27, 55, 9}, 
{97, 56, 165, 27, 55, 9}, 

}; 

static short OptSizeLUT [256 ] = { 



120, 

118, 

117, 

115, 

114, 

112, 

111, 

109, 

108, 

106, 

105, 

104, 

102, 

101, 

100, 

99, 

97, 

96, 

95, 

94, 

93 , 

92 , 

91, 



86 
85 
84 
84 
83 
82 
82 
81 
81 
80 
80 
79 
79 
78 
78 
77 
77 
77 
76 
76 
76 
75 



90 
90 
89 
88 
87 
86 



5 r s 
- 



75 
74 
74 
74 
74 
73 
73 
73 
72 
72 
72 
72 
71 
71 
71 
71 
70 
70 
70 
70 
70 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 



69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 
69, 



69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
69 
68 
68 
68 
68 
68 
67 
67 
67 
66 
66 
66 
65 
65 
65 
64 
64 
63 
63 
62 
62 
61 
61 
60 
60 
59 
59 
58 
58 
57 
57 
56 
56 
55 
55 
54 
54 
53 
53 
52 
52 
52 
51 
51 
50 
50 




49, 
49, 
48, 
48, 
48, 
47, 
47, 
46, 
46, 
46, 
45, 
45, 
44, 
44, 
44, 
43, 
43, 
43, 
42, 
42, 
42, 
42, 
41, 
41, 
41, 
40, 
40, 
40, 
39, 
39, 
39, 
38, 
38, 
38, 

}; 

static short OptDensi tyLUT [2 56] 

128, 

127, 

126, 

125, 

124, 

123, 

122, 

121, 

120, 

119, 

119, 

118, 

117, 

117, 

116, 

115, 

115, 

114, 

114, 

113, 



113 
112 
112 
111 
111 
110 
110 
110 
109 
109 
109 
108 
108 
108 
107 
107 
107 
107 
107 
106 
106 
106 
106 
106 
105 
105 
105 
105 
105 
105 
105 
104 
104 
104 
104 
104 
104 
104 
104 
104 
103 
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